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We propose a new approach to the Casimir effect based on classical ray optics. We 
define and compute the contribution of classical optical paths to the Casimir force 
between rigid bodies. Our approach improves upon the proximity force approxi- 
mation. It can be generalized easily to arbitrary geometries, different boundary 
conditions, to the computation of Casimir energy densities and to many other sit- 
uations. 
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1. Introduction 

The last 10 years have witnessed quite a revolution in the experimental 
techniques used to prove Casimir effect. 1 Casimir's original prediction for 
the force between grounded conducting plates due to modifications of the 
zero point energy of the electromagnetic field has already been verified to 
an accuracy of a few percent. Progress has been slower on the theoretical 
side. Beyond Casimir's original study of parallel plates, 2 we are only aware 
of useful calculations for a corrugated plate 3 and for a sphere and a plate. 4 
Simple and experimentally interesting geometries like two spheres, a finite 
inclined plane opposite an infinite plane, and a pencil point and a plane, 
remain elusive. The Proximity Force Approximation (PFA) was shown by 
Gies et al.^ to deviate significantly from their precise numerical result for 
the sphere and plane. Thus at present neither exact results nor reliable 
approximations are available for generic geometries. It was in this context 
that we recently proposed a new approach to Casimir effects based on 
classical optics. 6 The basic idea is extremely simple: first the Casimir energy 
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is recast as a trace of the Green's function; then the Green's function is 
approximated by the sum over contributions from optical paths labelled 
by the number of (specular) reflections from the conducting surfaces. The 
integral over the wave numbers of zero point fluctuations can be performed 
analytically, leaving a formula which depends only on the properties of 
the paths between the surfaces. This approach will give an approximation 
(though a surprisingly accurate one) which is valid when the natural scales 
of diffraction are large compared to the scales that measure the strength 
of the Casimir force. In practice this will typically be measured by the 
ratio of the separation between the conductors, a, to their curvature, R. 
It generalizes naturally to the study of Casimir thermodynamics, energy 
and pressure, to various b.c, to fermions, and to compact and/or curved 
manifolds. 

In the optical approximation the cutoff dependent terms of the Casimir 
energy can easily be isolated and shown to be independent of the separation 
between conductors. They therefore do not contribute to forces and can be 
dropped. 



2. Derivation 

Most studies of Casimir energies do not consider approximations. Instead 
they focus on ways to regulate and compute the sum over modes, ^ \^ w - 7 
These methods have proved very difficult to apply to geometries other than 
parallel plates. The main reason for this impasse lies in the requirement 
of an analytic knowledge of the spectrum of the Laplace operator for the 
given geometry. However, the knowledge of this spectrum for a family of 
boundaries (like a sphere facing a plane) would have the most non-trivial 
implications for the same family of quantum billiards and hence of classical 
billiards. 30 years of work on the ergodicity of classical billiards and their 
quantum counterparts suggest this task is hopeless. 8 

A numerical knowledge of the spectrum does not represent a reliable 
solution to the problem either. The force, indeed, is given by the small 
oscillatory ripple in the density of state numerically shadowed by the 'bulk' 
contributions which give rise to distance-independent divergencies. 

So we focused our attention on ways to get approximate solutions of the 
Laplacc-Dirichlct problem which are apt to capture the oscillatory contri- 
butions in the density of states, providing physical insights and accurate 
numerical estimates. 
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2.1. The Optical Approximation for the propagator 

There is a strong connection between the Casimir energy for a field <f> obey- 
ing Dirichlct boundary conditions (b.c.) on the boundary of the domain 
D and the propagator G(x',x,k) of the Helmholtz equation on the same 
domain with the same b.c. . The knowledge of the latter allows one to cal- 
culate the density of states p(k) and from this we can obtain the Casimir en- 
ergy by quadratures. Indeed, from the well-known definition of the Casimir 
energy in terms of a space and wave-number dependent density of states, 9 
p(x,k), 



where u(k) = C\Jk 2 + p 2 , and the density of states p(x, k) is related to the 
propagator G(x', x, k) by 



where the usual density of states is p(k) — J d N xp(x, k). 

We must choose G to be analytic in the upper-half A: 2 -plane; in the time 
domain (see later) this means we are taking the retarded propagator. 

The Casimir energy depends on the b.c. obeyed by the field <j> and on 
the arrangement of the boundaries, S = 3D (not necessarily finite), of the 
domain D. From the outset we recognize that E must be regulated, and will 
in general be cutoff dependent. We will not denote the cutoff dependence 
explicitly except when necessary, p and G are the familiar density of states 
and propagator associated with the problem 



(A + k 2 )^{x) =0 for x e D, 4)(x) =0 for x e S. (2.3) 



We can regard this problem as the study of a quantum mechanical free 
particle with K = 1, mass m = 1/2, and energy E — k 2 , living in the 
domain D with Dirichlet b.c. on dD. Dirichlet b.c. are an idealization 
for the interactions which prevent the quantum particle from penetrating 
beyond the surfaces S. This is adequate for low energies but fails for the 
divergent, i.e. cutoff dependent, contributions to the Casimir energy 10 
However, the divergences can be simply disposed of in the optical approach, 
and the physically measurable contributions to Casimir effect are dominated 
by k <~ 1/a, where a, a typical plate separation, satisfies \ ja <C A whit A 
being the momentum cutoff characterizing the material. So the boundary 
conditions idealization is quite adequate for our purposes. Following this 




(2.1) 



p(x,k) = — Im G(x,x, k). 



(2.2) 
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quantum mechanics analogy we introduce a fictitious time, t, and consider 
the functional integral representation of the propagator. 11 The space-time 
propagator G(x',x,t) obeys the free Schrodinger equation in D bounded 
by S. It can be written as a functional integral over paths from x' to x 
with action S(x',x,t) — j J Q dtx 2 . The optical approximation is obtained 
by taking the stationary phase approximation of the propagator G in the 
fictitious time domain 

G opt (x',x,t) = J2 K r(x',x,t)e iSAx ' x ' t) . (2.4) 

r 

The classical action is 

S r (x',x,t) = ^ X ; X ^ (2.5) 



At 



and K r is the van Vleck determinant 



With some manipulations 6 we can turn this determinant into 

where ./V is the number of spatial dimensions. This approximation is exact 
to the extent one can assume the classical action of the path S r to be 
quadratic in x',x. This is the case for flat and infinite plates. Thus the 
non-quadratic part of the classical action comes from the curvature or the 
finite extent of the boundaries, which we parameterize generically by R, 
d 3 S/dx 3 ~ 1/Rt. 

In fc-spacc the corrections hence will be O (1/kR), and the important 
values of k for the Casimir energy are of order 1/a, where a is a measure 
of the separation between the surfaces. Thus the figure of merit for the 
optical approximation is a/R. Certainly some of the curvature effects are 
captured by the van Vleck determinant but at the moment there is no 
good way to estimate the order in a/R of the corrections to the optical 
approximation (possibly fractional, plus exponentially small terms). This 
is topic for further investigation. 

Putting all together we find the space-time form of the optical propa- 
gator to be 

- £ (tf-^) (2.8) 
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When dealing with infinite, parallel, flat plates this approximation becomes 
exact. For a single infinite plate, for example, the length-squared of the only- 
two paths going from x to x' are l\ iroct = 1 1 x' - x \ \ 2 , Reflection = 1 1 x ' ~ * 1 1 2 » 
where x is the image of x. Both are quadratic functions of the points x, x' 
and the optical approximation is indeed exact. G opt (x' , x, k) is obtained by 
Fourier transformation and can be expressed in terms of Hankel functions, 
giving us the final form for our approximation 

G opt (x>, x,k)=^ (e?- 1 K) 1 ' 2 (I) Hf_ x (ki r ) , 

= ^G r (x',x,k), (2.9) 

r 

where A r is the enlargement factor (see 6 for details) 

M/,x, = ^ (2 . 10) 

and we have suppressed the arguments x and x' on l r and A r in (2.9). This 
can be thought of as a particular case of a more general result. 12 

2.2. The Optical Casimir energy 

The substitution of (2.9) into (2.2) and then in (2.1) gives rise to a series 
expansion of the Casimir energy associated with classical closed (but not 
necessarily periodic) paths 

E o P t= ]T Er, (2.11) 

paths r 

where each term of this series will be in the form of 

1 f°° 2k f 

E r = -hlm dkco(k)— d N x G r (x,x,k). (2.12) 

2 JO ^ JD r 

Here the integration has been restricted to the domain D r C D where the 
given classical path r exists. If the length of the path is bounded from 
below (this is the case for more than one reflection) the x and k integral 
can be switched safely. The fc-integral can be performed exactly for any N 
and /z, but it is particularly simple for the massless case, uj(k) = ck, 

( _ ir+ i r (N+l\ f N Al /2 

E r = nC 27r N/2 + l/2 r [-2- ) d X ^WT2- ( 2 - 13 ) 

This is the central result of our work and associates a Casimir energy con- 
tribution to each optical path r > 1. The series (2.11) has a very fast 
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convergence, usually 98% of the contribution is contained in the first 4 
terms, as we will see in the examples below. 

2.3. Divergencies 

The paths r = IB that bounce only once on a given body B must be 
treated with particular care because their contribution is divergent. The x 
and k integrals cannot be inverted without regulating the divergencies. To 
do so we insert a simple exponential cutoff in k. For a massless field the 
/c-integration in E\b can be performed* 1 giving 

ElB ~ "45? J DlB d x AlB (i + (W) 2 ) 3 • (2 - 14) 

Notice that for £isA > lwe reobtain the standard result, eq. (2.13). Hence 
it is convenient to rewrite this integral as 

Jo 1B ) 1B (1 + (W) 2 ) 3 ' 

(2.15) 

where Dib is the domain where the path does not exist. The first integral 
does not depend on the position of the other bodies and in the second we 
can take the limit A — > oo safely because £ib(x) > 2a > 0, where a is 
the minimum distance between the bodies. Hence we will write E\b = 
EiB,dw + EiB,fin, the finite part being the integral over Dib (notice the 
extra minus sign). 

The first term contains all the divergencies arising for 1\b — » 0. No- 
tice that since when £isA = \/3 the sign of the integrand changes, the 
divergence is negative rather than positive, as one could have argued from 
(2. 13). 15 Of course the bulk contribution to the vacuum fluctuation energy 
comes from the zero-reflection term, which is positive. 

Using the expression of Aib for a general surface with principal radii 
of curvature R a ,b we find 

Sis.div ~ -t^cA 3 - A 2 -^ftc / dS f i- + -M + (In A) . (2.16) 
»7T 327T J \R a R b J 

These terms do not contribute to the forces between rigid objects. The 
form of eq. (2.16) invites comparison with the work of Balian and Bloch. 13 
One finds agreement in the surface terms but not in the curvature terms. 



a For simplicity we specialize to N = 3 although the analysis is completely general. 
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3. Parallel Plates 

Parallel plates provide a simple, pedagogical example which has many fea- 
tures — fast convergence, trivial isolation of divergences, dominance of the 
even reflections — that occur in all the geometries we analyzed. We assume 
for simplicity that the two plates have the same area S. The relevant paths 
are shown in Fig. 1 where the points x and x 1 , which should be equal, are 
separated for ease of viewing. For the even paths lin(z) = 2na, n = 1, 2, . . ., 
independent of z (here z is the distance from the lower surface). For the 
odd paths £2n-i,a(z) = 2(n — l)a + 2£, where ( — z, a — z respectively if 
a = down, up and n = 1,2, . .. . For planar boundaries the enlargement 
factor is given by A n = 




2 4 6 2' 5 

1 3 

(a) (b) 

Figure 1. a) Even and b) odd optical paths for parallel plates. Initial and final points 
have been separated for visibility. 



The sum over even reflections, 

is trivial because it is independent of z. The result is the usual Dirichlet 
Casimir energy. 1 The sum over odd reflections, after being regulated by 
point splitting, before removing the infinite part, gives 

he r ,„^2, f a , 1 he 2nS 



E odd = — dSy / dz ? = -—^—5-- (3.18) 

2n J h>Jo (e 2 + (2z + 2na) 2 ) 2 16^ 2 e 3 V ' 

The divergence as e — > is precisely what is expected on the basis of the 
general analysis of the density of states in domains with boundaries. 13,14 
Moreover, since it is independent of a, it does not give rise to a force. 
From Fig. 1 b) it is evident why this must happen. The total sum of the 
off reflection contributions is just the integral of the one reflection path 
extended up to oo and hence docs not depend on a. 
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The fact that the odd reflections sum up to a divergent constant is 
universal for geometries with planar boundaries, and to a good approxi- 
mation is also valid for curved boundaries. Note also that the sum over 
n in cq. (3.17), converges rapidly: 92% of the effect comes from the first 
term (the two reflection path) and > 98% comes from the two and four 
reflection paths. This rapid convergence persists for all the geometries we 
have analyzed due to the rapid increase in the length of the paths. 

4. Sphere and Plane 

We calculated the Casimir energy for the sphere and the plane up to four 
reflections. E\ and E 3 can be found analytically while E 2 and E4 must 
be computed numerically. Henceforth a will be the distance between the 
sphere and the plate and R the radius of the sphere. 

Comparison with the parallel plate case as a — > suggests the error 
due to neglecting the fifth and higher reflections to be <~ 2%. Hence we 
have plotted our results as a band 2% in width in Fig. 2. Since the frac- 
tional contribution of higher reflections decreases with a, we believe this is 
a conservative estimate for larger a. 

The proximity force approximation has been the standard tool for es- 
timating the effects of departure from planar geometry for Casimir effects 
for many years 7 . In this approach one views the sphere and the plate as a 
superposition of infinitesimal parallel plates (i.e. the sphere is substituted 
by a stairlike surface). The resulting expression is 

plate _ ir 3 hcR 1 

^PFA --l U0a 2 1 + a/R - ( 4 - 19 ) 

It is custom to factor out the most divergent term of the Casimir force in 
the limit a/R — > as predicted by the PFA, in this case — 7r 3 ftci?/720a 3 , so 
to write in general 

F = _ f{ a AR 

• M iT 720a 3 ' 1 ' 

which is the definition of /. Modern experiments are approaching accuracies 
where the deviations of f(a/R) from unity are important. PFA predicts 
that 

/ PFA (a/ J R) = l-i| + o(^), (4.21) 
while the optical approximation data predict 

/optical (a/R) = 1 + 0.05 a/R + . (4.22) 
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Beyond the limit of small a/ R, one must notice that the optical ap- 
proximation to the Casimir energy and the data of Ref. 4 both fall like 
1/a 2 at large a/R. In fact both are roughly proportional to 1/a 2 for all 
a. In contrast the PFA estimates of the energy falls like 1/a 3 at large a 
and departs from the Gies et al. 4 data at relatively small a/R. For pur- 
poses of display wc therefore scale the estimates of the energy by the factor 
-U40a 2 /ir 3 Rhc. The results are shown in Fig. 2. The dominant contri- 




ve 1/64 1/16 1/4 1/2 12 4 

a/R 

Figure 2. Sphere facing a plane case. Comparison between different methods. Results 
for Ref. 4 (stars with error bars), data from 6, superseded by this work (triangles), optical 
approximation (thick grey line), PFA (broken line). 



bution, always greater than 92%, comes from the second reflection. The 
fourth reflection contributes about 6% for a/R <C 1 and less as a/R in- 
creases. The contributions of the first and third reflections are very small 
for all a/R. A relevant result, confirmed by the analytical analysis on the 
energy momentum tensor (within the optical approximation 16 ) is that the 
asymptotic behavior of E as a/R 3> 1 predicted by the optical approxi- 
mation is oc 1/a 2 . This is in contrast with the Casimir-Polder law which 
predicts E oc 1/a 4 at large a, the discrepancy is to be attributed to settling 
in of diffraction effects. 



5. Conclusions 

We have proposed a new method for calculating approximately Casimir en- 
ergies between conductors in generic geometries. Wc use an approximation 
imported from studies of wave optics that we have therefore named the 
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"optical approximation" . In this paper, we have outlined the derivation 
and applied it to two examples: the canonical example of parallel plates 
and the experimentally relevant situation of a sphere facing a plane. Our 
results are in agreement with the Proximity Force Approximation only to 
leading order in the small distances expansion. The first order correction 
is found to be different. This is of particular importance in the example 
of the sphere and the plane because the first order correction in a/R (a is 
the distance sphere-plate and R is the radius of the sphere) will soon be 
measured by new precision experiments. 

We would like to thank the organizers of the conference 'Continuous 
Advances in QCD 2004'. This work is supported in part by the U.S. De- 
partment of Energy (D.O.E.) under cooperative research agreement #DF- 
FC02-94ER40818 and in part by the INFN-MIT 'Bruno Rossi' fellowship. 
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